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Abstract 

We show how a precise knowledge of parton distribution functions, in particular those of 
the u and d quarks, can be used to constrain a certain class of New Physics models in which 
new heavy charged resonances couple to quarks and leptons. We illustrate the method by 
considering a left-right symmetric model with a W 1 from a SU(2)r gauge sector produced 
in quark- antiquark annihilation and decaying into a charged lepton and a heavy Majorana 
neutrino. We discuss a number of quark and lepton mixing scenarios, and simulate both 
signals and backgrounds in order to determine the size of the expected charge asymmetry. 
We show that various quark-W' mixing scenarios can indeed be constrained by charge 
asymmetry measurements at the LHC, particularly at \fs = 14 LeV. 

1 Introduction 

Many theories beyond the Standard Model (SM) predict existence of new particles that can be 
discovered at the LHC. Once they are discovered and basic properties such as masses measured, 
a next important step would be to determine the couplings of these new particles. Determination 
of these couplings will be crucial in unravelling the structure of the underlying physics models. 
While this is a non-trivial task, strategies to measure various important model-specific couplings 
have been proposed. 

In this paper, we discuss how precise knowledge of parton distribution functions (PDFs) 
could provide important information on a particular class of couplings. These are the couplings 
between quarks and charged resonances that might be produced at the LHC, i.e. via qq' — > X . 

The fact that the LHC is a proton-proton collider means that for many processes, positively 
charged (leptonic) final states will be more copiously produced than their negatively charged 
counterparts. Depending on the typical scales involved, the charge asymmetry ratio 



where N + (N_) are number of events with positive(negative) charge assignment, for a given event 
topology can be determined if the underlying hard process is knownQ In Ref. pj the variable R± 
was suggested as a diagnostic tool for the presence of New Physics (NP) with a SM background 
dominated by W + n jet processes^! An important advantage of using R± is that theoretical 
uncertainties tend to cancel in the ratio, leading to relatively robust predictions. 

Another feature of R± is that for (sub) pro cesses produced at the same scale, contributions 
from different incoming partons will contribute to the ratio differently, due to the non-universal 
flavour content inherent in the proton. For example, processes initiated by valence quarks u 
and d will tend to give R± > 1, since u(x) > d(x), while those initiated by the sea quarks c, s 
and b will lead to R± ~ 1. A measured value of R± can therefore lead to a constraint on the 
relative contribution from these initial state partons. Moreover, as the ratio is independent of 
the absolute normalisation of the couplings involved, it provides a genuine test of the structure 
of the parton coupling to the final state particle (s)q 

These ideas can be applied to NP models where new charged particles are produced at reso- 
nance at the LHC. The presence of a resonance implies that, in the absence of SM background, 
R± should track closely the ratio of 'weighted' parton luminosities 



ab 



dM v 



dy\V ab \ 2 f a {x u M v )f b {x 2 ,M v ), (2) 



where x\ = ^e y , x 2 = y , My = y/xix 2 s is the mass of the resonance, V ab is its coupling 
to partons a and b, and [y m i n , y m ax\ is the rapidity interval in which the resonance is produced. 
We denote such parton luminosity ratio by R±. At higher orders, additional QCD radiation 
will skew the parton luminosity prediction slightly, as will cuts on the decay products of the 
resonance. Nevertheless it should still be possible to obtain quantitative estimates of R± using 
appropriate Monte Carlo programmes. 

Charged resonances that might be discovered at the LHC include, for example, vector bosons 
in technicolor or extended gauge theories 0, Q , and charged sleptons in standard supersymmetric 
models with R-parity violation p|. In the same spirit, SM Higgs-stralung off a (off-shell) W boson 
is another interesting example, since we would expect a(W + H) > <j(W~H). 

For concreteness, we will focus in this study on a left-right (LR) symmetric model jsj with a 
heavy W from a SU(2)r gauge sector. The primary leptonic decay mode of the W is W — > INi, 
where iVj is a heavy Majorana neutrino, followed by three-body decay of Ni into a charged lepton 
and 2 jets. The Majorana nature of iVj leads to 'same-sign di-lepton (SSDL) plus 2 jets' events. 
This mode has not so far been searched for in the Tevatron (jpp) experiments (however see the 



CDF |9| and DO [lOj W searches in other decay channels). The discovery potential of this 
model at the 14 TeV LHC has been studied in Q We will see that the SM background 
can be strongly suppressed by applying cuts similar to those in the above studies, together with 



1 The paradigm process is of course Standard Model W ± production, for which R± ~ 1.33(1.42) at 14(7) TeV 
LHC in NLO QCD perturbation theory [3, 0]. Here the process is initiated primarily by ud and du scattering. 
2 For state-of-the-art calculations of W + n jet cross sections in the SM, the reader is referred to Refs. 0, 0j- 
3 Of course the statistics and hence the accuracy in determing R± will still be affected by the normalisation 
of the couplings involved. 
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the requirement that the di-lepton pair has the same charge. As a result, the experimentally 
determined value of R± will be correlated with the parton luminosity predictions. Note that the 
charge asymmetry ratio is only weakly dependent on the decay modes of the charged resonance. 
Therefore our discussion can be readily extended to other W models, for example the extended 
gauge theory model considered in j^l, 0] (the discovery potential of this model at the LHC is 
discussed in in addition to the models mentioned above. 

There exist models with additional constraining flavour structures, for example discrete left- 
right symmetries relating the SU{2) L and SU{2) R sectors. Without a concrete high-scale model 



implementation, the flavour/generation dependent couplings may be considered arbitrary [13 
Phenomeno logically, if W is relatively light, the quark mixing structure can be strongly con- 
strained by low energy observables such as flavour changing neutral currents (FCNC), for exam- 
ple from K — K° mixing [14J . We will illustrate to what extent different flavour structures can 
be distinguished using R±. 

The paper is organised as follows. In the next section we specify our W model and discuss 
the behaviour of the parton luminosity ratios R± given a set of representative quark flavour 
mixing matrices. In Section [3] we then perform a numerical study to show the relation between 
R± and R±, followed by a discussion on the extent to which R± can be used to distinguish quark 
flavour structures at 14 TeV LHC. Finally, we briefly comment on the prospect of using this 
observable in the early stages of LHC running. 

2 The W model 

The Lagrangian describing the coupling of the W to fermions (quarks and leptons) is given byQ 



C w > = -^%u i rP R d 3 W' u + + -Lu^PrNjW'- +h.c. (3) 

In the above, we have assumed that the W couples only to right-handed fermions with an 
overall coupling strength normalised to g. The quark and lepton mixings are given by V and 
U respectively. As our focus is on the effect of the structure of V on R±, we will assume, for 
simplicity, that g is the same as the SM electroweak coupling g and that U is diagonal, i.e. that 
the two charged leptons from W decay have the same flavour. This allows us to focus on same- 
sign, same-flavour leptonic final states. However for more general U, same-sign, different-Hecvouv 
final states can also be considered. 

In some (LR) symmetric models, V is related to the SM CKM matrix, Vckm- in the limit 
of (pseudo-)manifest LR symmetry, V = Vckm and the mass of W is constrained to be above 
a few TeV [lfjj], primarily by AF = 2 processes such as K° — K° mixing. However this lower 



bound can be relaxed by considering more general V structures, see for example Ref. [13]. A 



A general interaction Lagrangian is given in 15[ . 
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V 


Diagonal mixing 

o o 

CKM mixing 
Off-diagonal mixing 
1st generation only 
2nd generation only 


Vi = l 

VcKM — VcKM 

Vn (see eq. ®) 
(Vist)ij = V35n5ji 

(V2nd)ij = \ 35i26j2 



Table 1: Representative quark mixing matrices V in the charged current Lagrangian with a W, 
see Eq. fl3J). The normalisation of these matrices is such that \ Vij\ 2 = 3. 



recent study 18] suggests that if the quark mixing takes the form 

|V| ~ 1 , (4) 




then right-handed current contributions could lead to sizeable CP violating effects in B s mixing 



as suggested by recent results [19| , while remaining consistent with stringent AF = 2 observables 
for relatively low effective scales. We will take this quark mixing structure, from now on denoted 
by Vii, as one of our 'reference' structures. This can be compared with other representative 
structures, for example CKM-like mixing (Vckm), diagonal mixing (Vi) and when the W couples 
only to a particular generation (Vi st ,V2nd)- These quark mixing patterns are listed in Table l2lFI We 
will compare to what extent R± can distinguish these different mixing patterns. For simplicity, 
we assume no mixing between the SM W and the W . The decay mode W — > WZ is forbidden 
as a result. 

As discussed in the Introduction, once V is specified, the ratio of parton luminosities defined 
in Eq. [2] can be determined. As long as the events selected to determine the charge asymmetry 
ratio come dominantly from the W signal, R± should track the corresponding R± closely. We 
again note that this correspondence should depend only weakly on the event topologies and 
acceptance cuts. 

In Fig. [I] we display R± for the various quark mixings shown in Table [2] as a function of 
M w . In these plots, MSTW08 NLO PDF sets are used, with 68%cl PDF uncertainties also 
included. Once NLO pQCD corrections are accounted for, the PDF uncertainties are expected 
to be the most important theoretical uncertainty on the ratio predictions. The factorisation 
scale is chosen to be equal to the mass of the resonance. As we can see,the ratios for Vi st , V\ 
and Vckm are all quite similar, as they are dominated by contributions from the first generation 
quarks. However Vn shows a clear deviation from Vi and also V^^a when compared with the PDF 
errors J§ 

We can also understand the general behaviour of R± for different quark mixing patterns. 
As My/i increases, successively higher parton momentum fractions x are probed. While the 



5 The matrices are normalised according to J2ij l^il 2 = 3- 

6 Note that early LHC measurements of the charge asymmetry in SM production will very likely reduce 
the PDF errors further. 
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u/d ratio increases monotonically with x, the d/u ratio first increases and then decreases as x 
increases in the relevant range. The fractional decrease of the latter overcomes the increase of 
the u/d ratio, leading to an overall slight decrease of the parton luminosity ratio for Vi st , V\ 
and Vckm in the region My/' = 3 — 4 TeV. On the other hand, the luminosity ratios for the 
second and third generation partons are much flatter functions of x, leading to an increasing R± 
primarily controlled by the u/d ratio for Vn, and a R± ratio close to 1 for X^nd- 
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M w , (GeV) 




1000 2000 3000 4000 
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Figure 1: Variation of parton luminosity ratios R± with Mw at 14 TeV LHC. The quark mixing 
matrices V are defined in Table HJ The rapidity y of W satisfies \y\ < 2.5. The ratio is 
computed using the MSTW08 NLO PDF set |2j with factorisation scale fip = M W /. The error 
bars correspond to 68%cl PDF errors. 



In the next section, we will perform event simulations to quantify the correlations between 
R± and R±. We comment at this point that both values depend on the PDF sets used. For 



example, our simulation is performed using Herwig++2.4.2 20], interfaced with the LHAPDF 



library |2jJ, with MRST LO** as the default PDF set. Compared to MRST LO**, we find 
significant differences for the predicted values of R± and hence R± when using MSTW08 NLO 
PDFs (interfaced via the LHAPDF library j2l[). This can be understood by considering the 
differences in individual PDFS between the recent MSTW08 and earlier MRST LO** sets. Of 
course a detailed simulation of W production would require a proper NLO calculation using an 
up-to-date NLO PDF set such as MSTW08. We expect that the bulk of the NLO corrections 
comes from the differences between the LO and NLO PDF sets (see, for example, Ref. [l|), and 
so for this reason we will use MSTW08 NLO as our preferred PDF choice. Using other PDF 
sets, for example CTEQ6.6(NLO) |22j], will also lead to slightly different central values and PDF 
uncertainties. The discrepancies between different PDF sets will hopefully be resolved in the 
future. 
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3 Numerical study 



In this section we discuss the basic properties of the signal and backgrounds calculated using 
Herwig++2 .4.2. A set of cuts is then proposed to extract the signal from the background, before 
the prospect of quantitatively measuring R± is discussed. 

As discussed in the Introduction, the signal is same-sign di-leptons in association with 2 
hadronic jets. We assume that the heavy Majorana neutrinos are lighter than the W. Clearly, 
an accurate measurement of R± will require much higher statistics than a W discovery. We will 
therefore investigate a low mass range Mw> = 1 — 2 TeV. For simplicity, the masses of the 3 
heavy Majorana neutrinos are assumed to be degenerate, and are given by = \M W i. 

To identify a signal event, the basic strategy is to look for a high px SSDL pair and a pair 
of high pt jets coming from W — > INi — > lljj. The invariant mass of these four objects will 
lead to a resonance peak for the W . One of the two lepton-jet-jet combinations will also give a 
resonance peak for iVj. 

The dominant background is expected to come from tt. Compared with the studies in 



111 . [12], the requirement of SSDL efficiently suppresses this background. This is because at the 
hadron level, one of the leptons must come from B decay (the leptons from the two W decays 
naturally have opposite sign), which in general does not satisfy a lepton isolation requirement. 
Other potential backgrounds are from Wbbj and WZ( / y*)jj production. The former process 
will again be suppressed by a lepton isolation requirement. However its production cross section 
is relatively large, and therefore it demands more attention. The latter WZ (7*) jj process can 
lead to SSDL + 2 jet events if the wrong-sign lepton from Z(j*) — > l + l~ is not identified, for 
example when it falls outside the central tracking region. 

Note that since the tt background is produced predominantly by gg scattering, we have 
R± ~ 1 for this process, i.e. any final state involving charged leptons is symmetric under change 
of sign, whereas R± from Wbbj and WZ(^y*)jj is generally > 1 since these arise from qq' and 
qg scatterings. 

In terms of other possible backgrounds, we expect that by requiring SSDL, the Z(j*) back- 



ground discussed in [11|, [12| will be suppressed to a negligible level. The W W jj background 



will also be insignificant due to its small production rate [23] , l24j . Therefore these processes 
will not be discussed any further. We also do not attempt to study detector effects. For this 
reason, the QCD multijet background is not considered. A discussion of this background and its 



relevance for W discovery, including a detailed detector simulation, can be found in Ref. [11 
Note that any contamination of the W — > lljj signal from QCD multijet production will tend 
to reduce R± since this background is expected to be charge symmetric. 

To simulate the W signal, we implement the W model into Herwig++. The latter is also 
used to generate both the signal and the ti background. We use MSTW08 NLO PDFs together 
with default values for other parameters. The WZ( / y*)jj and Wbbj processes are simulated 



using ALPGEN[25], with Herwig6.510 26] then generating the parton showers. Parton level 
cuts appropriate to the experimental cuts discussed below are imposed to speed up the event 
generation, and we use 'out of the box' SM parameters for these processes. 
An event is selected for further analysis based on the following criteria: 
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• the leptons (jets) should lie within pseudorapitidy range \rj\ < 2.5 (4.5); 

• a lepton is isolated if E\ so < E^g = 5 GeV, where £^ so is the transverse energy in a cone 
of R = 0.4 surrounding the lepton; 

• the 2 hightest px leptons must have the same sign and flavour; 

• the 2 highest pt jets and the 2 highest pr leptons must be separated by AR > 0.4. 

In the above, an isolated lepton is assumed to be identified with 100% efficiency when p l T > 
10 GeV. Muons are assumed to be invisible when 1 77^ | > 2.5. Our jet reconstructions are 



performed using Fast Jet-2 . 4 . 2 |27|] using the anti-fc^ algorithm [28| with R = 0.4. 



We then apply further cuts to improve the signal/background ratio: 

1. select events with > 2 isolated leptons with p l T > 75 GeV and > 2 jets with pP T > 50 GeV; 

2. veto an event if a third lepton with p l T > 10 GeV is present which has the same flavour 
but opposite charge to the two hardest p l T leptons; 

3. retain an event when the invariant mass of the 2 leptons and 2 jets, mujj, satisfies 

0.7M^< m Ujj <1.2M^. 
Here Mffi is the value of the reconstructed W mass. 

Of course the backgrounds considered will have missing transverse energy If? already at the 
parton level. Furthermore leptons from B decays can have displaced vertices. These can be used 
as additional handles to further suppress the backgrounds. However we refrain from doing so here 
in an attempt to remain conservative. The leading-order multiparticle final-state backgrounds 
are known to suffer from scale uncertainties, and so the values quoted below should be treated 
as rough estimates only. 

Due to tight lepton isolation and hard lepton pt cut, we find that the Wbbj background 
is suppressed to negligible level. We provide a numerically more conservative estimate using 
MCFM [3] (at leading order). This parton level simulation does not model semi-leptonic decay 
of the 6s, rather a 'rule of thumb' efficiency factor eb ~ 1/200 is used for the probability of 
obtaining an isolated lepton of a given flavour (i.e. e or /x) from B decay suggested recently in 



29 



In Table [2] we display the effect of the selection cuts for M\yi = 1.0, 1.5,2.0 TeV with the 
quark mixing structures Vj and Vn- ^1,2,3 refer to the cross sections after applying the cuts 
1 (hard leptons and jets), 2 (opposite-sign, same-flavour veto (OSSFV)) and 3 (11 jj invariant 

mass constraint) defined above, including the ++ and contributions from both electrons and 

muons. In general, the signal cross sections for Vn are smaller than those for Vj, reflecting the 
different PDF combinations entering the processes via off-diagonal and diagonal flavour couplings 
respectively. Applying the M W / invariant mass cut further suppresses the backgrounds by a 
factor of a few, while retaining most of the signal as expected. Our more conservative estimate 
of the Wbbj background is still relatively small compared with the ti background. For the 
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Table 2: Cross sections (in fb) for different signal and background processes at 14 TeV LHC All 
4 same-sign, same-flavour di-lepton channels are included. &i,2,3 correspond to the cross sections 
after applying cuts 1, 2, 3 described in the text. e# denotes the efficiency in extracting an isolated 
lepton from a b quark and is estimated to be ~ 1/200 j29j. See Table [3] for the simulated R± 
values. 



WZ (7*) jj background, the OSSFV is shown to be efficient. In the following we will therefore 
neglect the Wbbj and W Z{^/*)jj contributions when performing the R± estimation. 

In Fig. [2] we show some sample kinematic distributions (normalised to 5 fb -1 ) from different 
processes. For clarity, only W' signals with V\ mixing and the dominant ti background are 
shown. We see clearly the W' and Ni resonance peaks from the invariant masses of lljj and 
Ijj respectively. The Ijj peak from the second hardest lepton is more pronounced due to our 
assumed relation = \M W i 1 which results in higher average pt for the lepton from the decay 
of the W' than that from the decay of the Ni. 

We also see that more pure samples can readily be obtained by varying the cut values, for 
example the px of the second hardest lepton, as a function of My/ 1 - Additional cuts, for example 
on the invariant mass of M^, may also be useful. Given that we already have fairly pure samples 
we refrain fro m op timising the cuts further. Other strategies for extracting the signal can be 



found in Refs. 11, 12 



As discussed earlier, the focus of this work is on the extraction of R±. In Tabled a comparison 
between R± and the simulated R± is shown. The cross-section ratios at the hard process level 
obtained from Herwig++2 .4.2 are well within the PDF uncertainties of the parton luminosity 
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Figure 2: Selected kinematic distributions for the same-sign dilepton events. The number of 
events are normalised to 5 fb~ . 
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71 /T 

M w > 


Vi 


V n 


tt wz{y)jj 


D TD 
±t± iX-|- 
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ii± JX± 


1.0 TeV 
1.5 TeV 
2.0 TeV 


2.12(4) 1.99(1) 
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Table 3: Comparison between parton luminosity ratio R± and observed charged asymmetry 
ratios R± (MSTW NLO) at 14 TeV LHC For the signal, Poisson uncertainties for the R±s at 
100 fb _1 and 1 a PDF uncertainties for R± are displayed in brackets. For the backgrounds, we 
show the statistical uncertainties based on the number of events generated. The Wbbj results 
are not shown due to too few statistics. 

ratio predictions, while the values of R± are systematically smaller than that of R±. The latter 
is due to the finite acceptance in lepton pseudorapidity and the fact that the average parton x 
for an up quark is higher than that of a down quark. The W' + , and hence the / + , will tend to 
have a broader rapidity distribution than that of the W'~ and l~, leading to a larger acceptance 
loss of the Z + . It is expected that while the extracted R± values will depend slightly on the cuts 
applied, they will nevertheless track R± closely. 

In Fig. [31 we compare R± and R± as a function of Myy, with quark mixings Vj and Vn as 
examples, at 14 TeV LHC with an assumed luminosity of 30 fb _1 . For simplicity, we have not 
included the background contributions to R±. From Table [3] we see that the latter have Re- 
values much closer to 1. Together with results in Table [2j including background contributions 
would lower the R± values slightly by 0(1 — 2)%. We see that for Mw> close to 1 TeV, the 
two mixing structures are clearly distinguishable. At this luminosity, the uncertainties on R± 
may be expected to come primarily from the PDFs. Statistical uncertainties become comparable 
at Myyi ~ 1.5 TeV. However at this luminosity, the ability to distinguish V\ from Vn becomes 
limited when My/' approaches 2 TeV. 

Finally, given the fact that in the My/ 1 = O(1.0) TeV region the W can be easily detected, we 
discuss the extent to which early LHC data (i.e. ~ lfb" 1 at 7 TeV) can measure R± . Measuring 
R± at different centre of mass energies will probe different parton x values, and hence could 
provide additional information to discriminate different PDF choices. 

In Table HI we show the results for the signal processes at 7 TeV, assuming an integrated 
luminosity of 1 fb -1 . Using the same cuts as the 14 TeV case, the backgrounds we have considered 
are again shown to be tiny (cr(tt) < lfb, a(WZ(j*)jj) < 0.2 fb and a(Wbbj) negligible after cut 
2). We see that while W could be discovered easily, the assumed total luminosity is only able 
to provide marginal discriminating power between V\ and Vn for Mw' below 1 TeV. However 
such low values of My/ 1 may have difficulty evading the bounds from EW precision observables 
already discussed. 
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Figure 3: Comparison between parton luminosity ratio R± and observed charged asymmetry 
ratios R± (MSTW NLO) at 14 TeV LHC as a function of Mw> ■ Results from two quark mixing 
patterns, V\ (upper, red curves) and Vn (lower, green curves) are shown. The colour bands 
correspond to 1 a PDF uncertainties, while the error bars are Poisson uncertainties on R± at 30 
fb" 1 . 



4 Conclusions 

In this paper we have demonstrated how by measuring the charge asymmetry ratio R± at the 
LHC, non-trivial information on the couplings of new heavy charged resonances to quarks may 
be obtained. We have focussed in particular on the W model, confirming expectations in the 
literature that in the M w > ~ O (TeV) region the physics backgrounds can be highly suppressed by 
applying appropriate cuts in the same-sign di-lepton + 2 jet channel. As a result, the measured 
R± value and the parton luminosity ratio R± are related in a directly correlated manner. Since 
the value of R± can be calculated precisely given a particular flavour mixing structure, we have 
established a method to constrain the quark mixing given a measured value of R±. Although 
we have considered overall event rates only, with higher luminosity it should be possible to 
study the rapidity dependence of the charge asymmetry ratio and this would provide additional 
constraints. 

Quantitatively, we have shown that at 14 TeV LHC three classes of representative flavour 
structures, namely V(i,cKM,ist), Vii and V2nd, can be clearly distinguished with realistic luminosi- 
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Table 4: Cross sections (in fb), R± and R± for different W models at 7 TeV LHC All 4 same- 
sign, same-flavour di-lepton channels are included. 03 corresponds to the cross section after 
applying cut 3 described in the text. The brackets for R± are ±1 a PDF uncertainties, whereas 
those for R± are Poisson uncertainties at 1 fb -1 luminosity. After applying cut 2, the background 
processes are suppressed below 1 fb, and hence their values are not shown. 



ties. However the prospect at 7 TeV LHC is less optimistic due to the lower overall luminosity 
as well as the larger PDF uncertainties in the higher parton x region being probed. 

A knowledge of quark flavour mixing structure will clearly have important implications on 
EW precision observables. Interestingly, if W is relevant for neutrinoless double beta decay 
(Ov/3/3), then measurement of i?-t0will be crucial in estimating its contribution to 0u(3(3, as only 
the first generation quarks participate in this low energy process. The same consideration will 
apply in relating an observation of resonant selectron production in R-parity violating SUSY 



models to its contribution to 0i/f3f3 30 



Finally, the measurement of R± can be extended to signals other than resonance production. 
We will study the application of this observable in other NP models, for example SUSY, in 
forthcoming publications. 
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